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MOVING TABLE MRI WITH SUBSAMPLING 



The invention relates to a magnetic resonance (MR) method for forming an 
image of an object from a plurality of signals sampled in a restricted homogeneity region of a 
main magnet field of the magnetic resonance imaging apparatus, whereas a patient disposed 
on a table is moved continuously through the bore of the main magnet, such that an image is 
5 formed over a region exceeding largely the restricted region of a main magnet field of the 
magnetic resonance imaging apparatus. 

The invention also relates to an MR apparatus and a computer program 
product for carrying out such a method. 

Since it is quite uncomfortable for a patient being disposed in a large bore of 
10 the main superconductive magnet of an MR imaging apparatus, there is a trend to use shorter 
magnets while moving the patient table or couch through the bore of the magnet. Such an 
apparatus is described e.g. in EP-A-1 024 371. In an MR imaging apparatus excitation pulses 
are applied to a restricted region of the magnet bore, in which the field is uniform. The data 
samples collected are Fourier transformed to form a volumetric image of the restricted 
15 region. A motor continuously moves a patient couch so that a region of interest passes 

through the region of good field. The collected data samples are corrected to compensate for 
the motion so that a volumetric image is formed of greater length than that of the restricted 
region. 

The MR imaging as mentioned above takes quite a long time since it is 
20 performed with conventional imaging sequences and even an overlap between successively 
imaged volumes is preferred in order to discard artefacts produced by errors in the central 
region. Obviously, the r.f. transmit coil and the r.f. receive coil are disposed in the reference 
frame of the main magnet, i.e. in its direct vicinitj'. However, imaging by body coils disposed 
on different places of a patient was not an issue in this reference. Anyway, it is known that 
25 using body coils may cause additional artefacts which cannot be suppressed in an easy 
maimer. 

Therefore, it is an object of the present invention to improve the MR imaging 
method with a relative short main magnet and a moving bed in such a manner, that imaging 
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will be performed in a much shorter time period. Another object of the present invention is to 
provide an MR apparatus and a computer program product for carrying out the method. 

These objects are achieved by a method as claimed in Claim 1, by an MR 
apparatus as claimed in Claim 1 1 and by a computer program product as claimed in Claim 
5 13. 

The invention has the main advantage, that shorter main magnets with a 
restricted homogeneity region can be used whereas data can be sampled in a much shorter 
time due to undersampling and therefore the scan time is largely reduced in relation to 
conventional MR imaging with a moving table or bed. Another advantage of the present 

10 invention is that the direction of the field-of-view can be oriented in any direction^ which 
means that e.g. flow artefacts can be oriented in an appropriate manner or to orient single 
slices to some preferred orientation with respect to the tissue structure of the patient. 

These and other advantages of the invention are disclosed in the dependent 
claims and in the following description in which an exemplified embodiment of the invention 

15 is described with respect to the accompanying drawings. Therein shows: 



Fig. 1 diagrammatically a magnetic resonance imaging system in which the 
present invention is used, 
20 Fig. 2 an example of a acquisition scheme, 

Fig. 3 an acquisition scheme in k-space with a moving bed, and 

Fig. 4 a simplified presentation for explanation of the present invention. 

Basic principles and notations 

25 In Fig. 1 a patient 1 is placed fixedly to a moving table or bed 2, A set of local 

coils LCi to LCs is located fixedly to the patient 1, thus the arrangement "patient V\ "bed 2" 
and "local coils LCi to LCs" are disposed in a fixed relationship to each other in a coordinate 
system 3 which is related to the patient. A vector x denotes the position relative to this 
coordinate system or fi-ame of reference 3. Further a main magnet 4 is provided with a bore 5, 

30 in which the bed 2 can be moved to and fi-o in direction of arrow 6. The magnet/gradient 
system 7 comprises the main magnet 4 and a number of receive coils 8 - generally one or 
none - and a transmitter coil (not depicted). These receive coils 8 and transmit coils are 
denoted as being "global". The elements of the magnet/gradient system 7 are all related in 
another fixed relationship which is defined by another coordinate system 9 related to the 
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main magnet 4. It is further assumed that the receive coils 8 and the transmit coils are large 
and have smooth sensitivity patterns in the direction of the bed motion (arrow 6). The magnet 
coordinate system 9 is offset by a distance A relative to the patient coordinate system 3. The 
property of the patient 1, e.g. the proton density of his body, is obviously given in the patient 
5 coordinate system 3, expressed as p(x) . Ideally, the same holds for the receive-sensitivity 
pattern S^(x) of one of the local coils LCi, where i = 1 to 5. On the other hand, the sensitivity 
pattern of the receive coil 7 is written as g/(x - A) . 

Now the expression "receive situation" is introduced, which means a block of 
measurements contiguous in time, during which the table 2 has travelled a relatively small 

10 distance (e.g. one half or one third or something in between) of the "homogeneity size" of the 
main magnet 4. During a single receive situation all magnetisation preparations (like 
presaturation slabs, volume selections etc.) "travel" in principle with the patient 1. Also 
during a single receive situation, all of the interesting k-space is covered, except for any 
(potentially significant) undersampling. The full scan over k-space covers K receive 

15 situations. In principle, the MR acquisition sequence is identical over all receive situations, 
except that the "magnetisation preparations" may change place and that there may be some 
offset in the acquired k-space grid. This is however a refinement of the concept. 

During any receive situation, every coil LCi to LCs receives information from 
the patient 1. And so (I -M) sets of information are delivered, where I is the number of local 

20 receive coils and J is the number of global receive coils. Some of these sets may be 

practically zero, e.g. if a receive coil LCi to LCs is completely outside the homogeneity 
volume of main magnet 4. There are K receive situations. Thus, in total, information has been 
gathered in K * (H- J) "receive instances". 

25 Basic idea of the invention 

The clinical doctor wants information from a relatively large volume of 
interest of the patient I. Data is acquired while the table or bed 2 is moving relative to the 
main magnet 4. The displacement A thereof relative to the magnet 4 is known at any time. 
The whole of the MR acquisition is split up in a number of receive situations. During each 

30 situation, the excitation and presaturation profiles "travel" with the patient 1 . The acquired 
data are phase-corrected for the offset of table 2 at which they are acquired. During a receive 
situation, information from the patient 1 is acquired using any MR imaging sequence of any 
orientation. In practice it may be useful either to orient artefacts (e.g. flow artefacts) in an 
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appropriate direction, or to orient single slices to some preferred orientation with respect to a 
specific human-tissue structure. However, the information may be seriously undersampled. 
An example of this method is sketched in Fig. 2, For every receive situation, we acquire, 
nominally the volume indicated by a box 11. This means that a multitude of points 12 will 
5 fold onto the results of this receive situation. This multitude may be a one-dimensional row, a 
two-dimensional row or even a continuum of folding locations which is the case in e.g. spiral 
acquisition. 

Now, every receive situation has the same subsampling pattern. The only 
difference between different receive situations is the weighing with respect to the properties 
10 of following selected factors: 

1. the excitation or saturation profile (with respect to the patient 1), 

2. the transmit pattern of the transmit coil, 

3. the receive patterns of the fixed receive coils 8, and 

4. the frequency-response pattern of the receiver system. 

15 Concerning items 2 and 3 it is remarked that in theory, with a moving table 2, 

these patterns varies also within a situation - but if the motion per situation is not exceedingly 

large, this may be neglected. 

In addition, every receive coil (whether local or global) has its own sensitivitj' 

pattern. So everj^ receive instance, which is the combination of a receive coil and a receive 
20 situation, "sees" a different "overall pattern", which is a product of the coil-sensitivity pattern 

and the combined pattern (excitation, saturation, transmit and frequency-response) of the 

receive situation. 

The multitude of folding points 12 is unfolded using all the "overall patterns" 
of every receive instance. This is done in a manner as in the SENSE method (see e.g. K. 

25 Pruessmann et. al. in Proc. ISMRM, 1998, abstracts pp. 579, 799, 803 and 2087) - although, 
formally, the unfolding matrix may be larger. Nevertheless, the system of equations may be 
reasonably stable: e.g. for each situation, many folding points 12 may practically "see" a zero 
pattern, because they fall outside the excitation slab. It is further to be noted that the table 
motion is flexible: it may be linear, non-linear and even non-monotonic or two-dimensional, 

30 The restriction is that it must not travel excessively during one receive situation, e.g. no more 
as half the homogeneity volume. 
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Measuring the patterns 
Calibration 

In SENSE, a map of coil sensitivity profiles must be available. Such a map can 
be provided in a reference or coarse calibration measurement with a well-defined phantom. 

5 Another possibility is using a low-resolution map of the ratio of signals from an array coil 
element and jfrom a reference coil of the patient on the table. In the present invention a 
similar reference measurement is provided. Although data can generally be acquired at any 
table-motion pattern, including two-dimensional table motion, here a linear one-dimensional 
movement will be used for data acquisition. 

10 For a better understanding of the involved principles, first a review: 

1 . At each position of the patient ( 3c ) some position-dependent tissue-properties 
(all data about proton density, Ti, T2, flow, diCEiision...) are known 

2. Also at each position of the patient ( 3c ) the sensitivity pattern (3c) of the 

local coils LCi can be obtained. Although there might be some dependence on A due to 
15 coupling to global coils, this is neglected here. 

3. The pattern of all the coils (global or local, transmit or receive) may be 
influenced by dielectric resonances, particularly at main magnet fields of 3T and higher. It is 
assumed that this effect can be seen as purely multiplicative on the coil sensitivity of all coils 
- an effect that is only dependent on the patient and the position therein, so not on the coil or 

20 coil-position. For a main magnet field of 1,5T the effect is negligible. The dielectric 

resonance effect is denoted as D(x) . Each coil actually "sees" with a sensitivity D(x) • S^ix) . 

4. The global coils also have their sensitivity pattern, which is approximated as 
being independent on the actual position of the patient within the magnet frame (except, of 
course, for dielectric resonance effect, which is explicitly accounted for). This sensitivity 

25 pattern is expressed as gj (x + A) . 

5. The total receive-pattem also may contain effects of receiver frequency 
response, denoted as F(x + A) at a given readout gradient direction. 

For a given measurement type, the transverse magnetisation is a non-linear 
effect depending on both the transmit pattern - i.e. on the position hi the magnet coordinate 
30 system 9, if a global coil is used - and the tissue propertj' - i.e. on the position in the patient 
coordinate system 3. Since this is very difficult to handle, an approximation is necessary, 
which can be achieved by using an FFE measurement with very small tip angles (about 5 
degrees or less). In that case, the transverse magnetisation can be approximated as 
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fix) • aix + A) , where f{x) is the "spin density" function, which is a tissue property for a 
given measurement (in the proposed case, it mainly depends on proton density), a is the tip 
angle, which is approximated as being proportional to D{x) • g,,ansmit • (5 + A) . 

5 Problems 

For the reconstruction of the proposed method as described in more detail 
below, ideally one would like to know all mformation of S^ix) , gy(3f + A) and Z)(x) (but 
the last one is not crucial). The value of F(x -i- A) is assumed to be known beforehand, and 
the value of f{x) is (at least in principle) not an interesting result of the calibration 
10 measurement. 

The problem is that not all these parameters are measured. For a given table 
position A, the following is measured: 

- the signal from local coils: F(x + A) • 5, i^Wx),g,^^^,^ (3c + A).Z)(5)./(jc) . 

- the signal from global coils: F(x + A) • g/^ + A).D(x).g,,^„,„„., (jc + A).D(3c)./(3c) . 

1 5 A plain division of those two signals results in (3c) / (3c + A) - which is of no direct use if 
gj(x + A) is not uniform. 

First solution 

The patterns of the global coils are assumed to be known ( i.e. at least one of 
20 them). In practice, that function may be load-dependent, but this can be overcome by the fact 
that there is a calibrated load-dependency of the sensitivity function, and the load can be 
easily measured. 

A number of discrete reference-scan segments are acquired. During such a 
segment, the table stands still. Each segment results in a full "image". The table is then 
25 stepped over a number of segments, so for a number of different values of A . The 

knowledge of gy(^ + A) allows calculating the matrix iS;(jc) . Mathematically, each segment 
would give the result over the full extent of 3c , but in practice, each segment will provide 
accurate results only over a sub-range of all 3c . These partially overlapping ranges can be 
combined using least-squares fits, resulting in the full map of S^ix) . 



30 
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Second solution 

The three-dimensional measurement is oriented in such a way that the 
frequency-encoding is parallel to the table motion direction - the z-direction. The sampling 
bandwidth is very high in that direction. The other two directions are phase-encodings, one of 
5 which is "fast" compared to table motion. The "slow" phase-encoding direction is called ky. 
In the patient coordinate system 9 the main magnet 4 is slowly moving with respect to the 
patient. During the motion of the magnet, profiles are acquhed with linearly increasing ky. In 
Fig. 3 a line segment 13 indicates the portion of the patient 1 that is detected by the 
homogeneity volume of the main magnet 4. 

10 L is the distance that the main magnet 4 covers for a full set of phase-encoding 

profiles. The table speed is arranged in such a way that L is maximally half of the 
homogeneity volume - but less (e.g., 1/4 ) is preferable. In Fig. 3 L is chosen at 1/2. The 
consequence is that every profile has been sampled for at least two different positions of the 
main magnet 4 (except for the outer edge regions, which are disregarded here). For example a 

15 profile ky has been measured for both the displacement Aj and for the displacement , This 
allows to approximate the profile that would have been measured if the main magnet 4 were 
at position (as long as this value is in between Ai and Aj ). This can be done by 

interpolation (e.g. a Imear interpolation if two profiles are available, or a cubic spline if there 
are four profiles). That can be done for every profile, which allows a complete reconstruction 
20 of the sensitivity-calibration measurement at magnet-offset A^^^ . 

However, this can also be done for any value of Are/ . That principle allows 

acquiring: 

- the signal from local coils: F{x + A) • (x) - D(x) • g^,„,,,, (x -f A) • D(x) • f(x) . 

- the signal &om global coils: F(x + A) • g (x + A),D(x),g,^^„^^^^, (x + A).D(x),f(x) . 

25 This holds, in contrary to the "problem section, not only for one value of A but for all values 
thereof. This allows to: 

- separate the signal from local coils into F(x + A) • gt,^,^^u + ^) and 
S^ix).D(x).D(x).f(x), and 

- separate the signal from global coils into F(x + A) • gj (x + A).g,^„^„„., (x + A) and 
30 Dix).Dix)J{x) . 
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This allows the derivation of S^x and of gj(x + A) . Actually, if F is known, then also 

Stransmu + ^) could be derlved, but unfortunately not the dielectric resonance D{x) - but 
that is not crucial for the sequel. 

5 Receiving and reconstructing the actual data 

In the section "basic idea of the invention" the receiving and reconstruction of 
data has been explained to a great extent. The data is acquired using relatively large steps 
between profiles, resultmg in a (nominally) small "folding volume", i.e., into lots of folding. 
For non-cartesian sequences this means that data is sampled relatively sparse. The full extent 
10 of k-space is acquired during one scan situation, and the displacement during one scan 
situation is a fi-action (e.g. 1/3) of the homogeneity volume. Every incommg sample is 
acquired at a very specific value of A^^,^^ . This is taken into account by multiplying every 

sample of the incoming data with exp( - /^JA^^^^^/^ ) . Here, T is the coordination 

transformation from patient coordinates to scan coordinates, which is relevant for oblique 
15 scanning. That operation "displaces" the acquired data to the centre of the patient, even if that 
is way out of the homogeneity volume. Obviously, one can simplify things by correcting with 
a fixed Aprqfi/e per profile. 

The scan situation has been acquired at an average offset A^,^ . This is relevant 

to estimate: 

20 - the transmit profile for that situation: (x) = (x + A^ J , 

- similarly for the global-coils receive profile, and 

- similarly for the frequency response function for that situation. 

The excitation profiles and presaturation profiles are fixed to the patient 
coordinate system, i.e. travelling with the table 2 as seen fi-om the main magnet 4. So these 

25 are known in x . Also the local coil profiles are known in x (see calibration section). All 
these coil profiles can be multiplied into an overall receive pattern for each coil. In total, if 
there are K receive situations, then information has been gathered hi K«(I + J) receive 
instances. Each instance has a different overall receive pattern. This allows for a SENSE-like 
reconstruction, if there are less than K '^(I + J) folded points over the entire patient 1. In 

30 principle, the acquisition stays the same over all receive situations. As a refinement, the 

phase encoding may be offset by a small amount. This will give some extra phase-encoding, 
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on top of all other mentioned encodings, which will improve the stability of the SENSE- 
reconstruction. 

Geometric correction 

5 All of above reasoning assumes a perfectly linear gradient system. It assumes 

that a table-displacement of A will cause (a) some different weighting, which has been 
extensively treated above, and (b) a phase modification of exp(/A:.rA,^,„^/^ ) to a sample. 
Unfortunately, geometric distortion is another significant complicating factor. An image 
acquired at offset Aj is geometrically distorted when compared to an image acquired at offset 

10 Aj , even if all weightings and the displacement ( Aj - A, ) have been accounted for. 

This problem can be solved since the distortions tend to be small where the 
excitation and transmit profiles tend to be large. This can be accomplished (a) by making 
sequences in such a way that there is not too much excitation in geometrically 
inhomogeneous areas, and (b) by designing the system with low transmit-sensitivity in these 

15 areas. The solution is an iterative reconstruction, the explanation thereof is given visually in 
Fig. 4. The main idea is, in essence, to first reconstruct the resulting image by disregarding 
the presence of geometric distortion. Of this result, all image regions are distorted in as many 
ways as there are scan situations. The distorted minus the undistorted versions provide an 
estimate of "folding distortion-errors". If the contributions (or weights) of every receive- 

20 instance to the final image are known, the estimate of the folding distortion-errors can be 
extracted from the result. 

In Fig. 4 (a) a very simplified patient 15 is depicted which diagrammatically 
consists of two parts at different positions: a first part shown as block 16 and a second part 
shown as heart 17. Fig. 4 (b) shows the imaging while the block 16 is largely in the centre of 

25 the main magnet and the heart 18 is more remote and thus distorted. As shown in Fig. 4 (c) 
the distorted heart 18 folds onto the block 16, due to undersampling. Yet, the heart 18 is 
fainter, since it has seen less excitation. In another receive-situation according to Fig. 4 (d), 
the heart 17 is in the centre of the main magnet 4 and the block 19 is more remote and thus 
distorted. Due to undersampling, this is measured in that receive-situation, as shown in Fig. 4 

30 (e). So the images in Fig. 4 (c) and Fig. 4 (e) are actually measured. Fig. 4 (f) shows the 
result after SENSE-unfolding, which consists of two linear combinations of Fig. 4 (c) and 
Fig. 4 (e). It is imperfect, e.g. the right part of it consists of the block 16 plus an unsubtracted 
heart 18 plus a "wrongly subtracted" heart 17 (dashed) and a second-order error (dotted). In 



wo 2004/077086 



10 



PCT/IB2004/050131 



Fig- 4 (g) is shown a purposely distortion in the same manner as in (b), under the condition 
that (f) is available. Fig. 4 (h) shows the final result where part of (f) is subtracted from its 
folded counterpart, and part of (g) is added, resulting in only second-order artefacts. 
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CLAIMS: 



L A magnetic resonance imaging method for forming an image of an object from 

a plurality of signals sampled in a restricted homogeneity region of a main magnet field of a 
magnetic resonance imaging apparatus, wherein a patient disposed on a table is moved 
continuously through the bore of the main magnet and spins in a predetermined area of the 
5 patient are excited by an excitation pulse from a transmitter antenna, such that an image is 
formed over a region exceeding largely the restricted region, characterized in that data is 
undersampled in the restricted region by means of at least one receiver antenna in a plurality 
of receive situations being defined as a block of measurements contiguous in time having 
preserved magnetisation and presaturation conditions within the excited area of the patient, 
10 and fold-over artefacts due to said undersampling are unfolded by means of the sensitivity 
pattern of the receiver antenna and/or the properties of selected factors determining said 
receive situations. 

2. A method as claimed in Claim 1, characterized in that the fold-over artefacts 
1 5 are unfolded by means of the excitation profile as selected factor. 

3. A method as clahned in Claim 1 or 2, characterized in that the fold-over 
artefacts are unfolded by means of the magnetisation and presaturation profile as selected 
factor. 

20 

4. A method as claimed in one of Claims 1 to 3, characterized in that the fold- 
over artefacts are unfolded by means of the frequency-response pattern of the receiver as 
selected factor. 

25 5. A method as claimed in one of Claims 1 to 4, characterized in that during 

sampling of data within the restricted region the table will be moved over at most one half of 
the size of the restricted region. 
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6. A method as claimed in one of Claims 1 to 5, characterized in that at least one 
global receiver antenna disposed in a fixed relationship to the main magnet system and a 
plurality of local receiver antennae disposed in a fixed relationship to the patient on the table 
are provided. 

5 

7. A method as claimed in one of claims 1 to 6, characterized that a reference 
scan is provided for obtaining the sensitivity pattern of the receiver antenna. 

8. A method as claimed in claim 7, characterized in that data is acquired in a 
10 number of discrete reference scan segments at different table positions in the restricted 

region, whereas during each reference scan the table stands still, and from the data obtained 
by the reference scans the sensitivity patterns of the local coils are calculated. 

9. A mefliod as clauned in claim 7 or 8, characterized in that the table is moved 
15 at a speed which is less than half of the homogeneity volume of the main magnet over the 

scan time of k-space, whereas k-space is continuously scanned in a row-to-row manner, and 
the sensitivity profiles of the local coils are calculated by interpolation of the measured 
profiles at different table positions. 

20 10. A method as claimed in any of claims 1 to 9, characterized in that the sampled 

data is reconstructed in an iterative manner, in that data sampled largely offset of the centre 
of the main magnet, which is folding-in on the data sampled in the centre of the main magnet, 
is purposedly distorted such that a undistorted image is reconstructed and subtracted from the 
fold-in image. 

25 

11. A magnetic resonance imaging apparatus for obtaining an MR image from a 

plurality of signals comprising: 

- a main magnet, 

a transmitter antenna for excitation of spins in a predetermined area of the patient, 
30 - at least one receiver antenna for sampling a plurality of signals in a restricted 
homogeneity region of the main magnet field, 

- a table for bearing a patient, 

- means for continuously moving the table through the bore of the main magnet. 
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- means for adapting the frequency profile of the transmitter antenna and the frequency 
profile of the receiver antenna to the position of the continuously moving table relative to 
a reference point, 

- means for sampling the data in the restricted region in a given receive situation being 
defined as a block of measurements contiguous in time with preserved magnetisation and 
presaturation conditions in the excited area of the patient, and 

- means for unfolding fold-over artefacts due to undersampling by means of the sensitivity 
pattern of the receiver antenna and/or the properties of selected factors determining said 
receive situations. 

12. Apparatus according to Claim 11, fiirther comprising at least one global 
receiver antenna disposed at a fixed relationship to the main magnet system and a plurality of 
local receiver antennae disposed in a fixed relationship to the patient on the table. 

13. A computer program product stored on a computer usable medium for forming 
an image by means of the magnetic resonance method, comprising a computer readable 
program means for causing the computer to control the execution: 

- creating a main magnetic field by a main magnet, 

- excitation of spins in a predetermined area of the patient by a transmitter antenna, 

- sampling a plurality of signals in a restricted homogeneity region of the main magnet 
field by at least one receiver antenna, 

- continuously movmg a table bearing a patient through the bore of the main magnet, 

- adapting the frequency profile of the transmitter antenna and the frequency profile of the 
receiver antenna to the position of the continuously moving table relative to a reference 
point, 

- undersampling the data m the restricted region in a given receive situation being defined 
as a block of measurements contiguous in time with preserved magnetisation and 
presaturation conditions in the excited area of the patient, and 

- unfolding fold-over artefacts due to said undersampling by means of the sensitivity 
pattern of the receiver antenna and/or the properties of selected factors determining said 
receive situations. 
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